NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CALIFORNIA 


THESIS 


EFFECTS  OF  MOSAICITY  ON  PARAMETRIC 
X-RADIATION  PRODUCTION 
by 

William  Morgan  Buckingham 
and 

Lisa  Renee  Ivey 
December,  1994 

ThgsisAdyisor _ Xavier  K.  Mamyama 

Approved  for  public  release;  distribution  is  unlimited. 


19950427  031 


REPORT  DOCUMENTATION  PAGE 


Fomi  Approved  0MB  No.  0704-0188 


Mc  reporting  burden  for  this  coflection  of  information  is  estimated  to  average  ]  flour  per  response,  including  the  time  for  reviewing  instruction,  searching  existing  data  sources,  gathering  and 
naintarning  the  data  rteeded,  and  completing  and  reviewing  the  coflection  of  information.  Send  comments  regarxling  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including 
uggesSons  for  reducing  this  burden,  to  Washington  Headquartcn  Services,  Directorate  for  Information  Operations  and  Reports,  121 5  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA 
2202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188)  Washington  DC  20503. 


1 .  AGENCY  USE  ONLY  (Leave  blank) 


2.  REPORT  DATE 

December,  1994 


».  TITLE  AND  SUBTITLE 

EFFECTS  OF  MOSAICITY  ON  PARAMETRIC  X-RADIATION 
PRODUCTION 


AUTHOR(S) 

Buckingham,  William  M.  and  Ivey,  Lisa  R. 


.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School 
Monterey  CA  93943 


.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


3.  REPORT  TYPE  AND  DATES  COVERED 

Master’s  Thesis 


5.  FUNDING  NUMBERS 


8.  PERFORMING  ORGANIZATION 


10.  SPONSORING/MONU  ORING 


1 1 .  SUPPLEMENTARY  NOTES  The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official  policy  or 
position  of  the  Department  of  Defense  or  the  U.S.  Government. 


I2a.  DISTRIBUTION/A VAILABILITY  STATEMENT  12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  is  unlimited. _ 

13.  ABSTRACT  (maximum  200  words) 

Previous  measurements  of  parametric  x-radiation  (PXR)  spectra  from  mosaic  graphite  observed  yields  that  greatly 
exceeded  theoretical  predictions.  At  the  Bragg  scattering  condition,  the  interaction  of  relativistically  charged  particles  and 
the  atomic  planes  of  a  crystal  produces  scattering  of  virtual  photons.  This  is  known  as  PXR.  The  discrepancies  between 
actual  yields  and  theoretical  predictions  might  possibly  have  been  caused  by  the  mosaic  spread  of  the  atomic  planes  of  the 
crystal.  Re-analysis  of  previous  data  has  shown  a  miscalculation  of  the  effective  target  thickness  used  in  earlier  work.  There 
was  also  an  error  in  the  calculation  of  the  theory  that  was  used  for  comparison.  Correction  of  the  two  errors  brought  the 
data  and  theory  in  closer  agreement,  but  the  mosaic  spread  does  change  the  intensity  of  the  PXR. 

In  this  thesis,  three  samples  of  thick  mosaic  carbon  graphite  crystal  with  differing  mosaic  spreads  are  compared  to 
determine  the  effect  of  the  mosaicity  on  PXR  yield.  By  obtaining  "rocking  curves"  where  the  target  orientation  angle  is 
changed  with  respect  to  the  incident  electron  beam  and  observation  angle,  it  was  found  that  the  mosaic  spread  docs  have  an 
effect  on  the  PXR  spectral  distribution.  This  effect  is  quite  significant  and  can  be  seen  when  viewing  the  "rocking  curves"  of 
the  samples.  As  the  mosaicity  is  increased,  the  angular  distribution  was  spread  out  and  the  coherence  condition  was  relaxed. 


14.SUBJECT  TERMS 

parametric  x-radiation,  rocking  curves,  mosaic  spread 


15. NUMBER  OF  PAGES 
_ 89 

16. PRICECODE 


I7.SECURITY  CLASSmCATION  18.SECURITY  CLASSIFICATION  1 9.SECURITY  CLASSIFICATION  20  LIMITATION  OF 
OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT  ABSTRACT 


Unclassified 


Unclassified 


Unclassified 


/ 


11 


Approved  for  public  release;  distribution  is  unlimited. 


EFFECTS  OF  MOSAICITY  ON  PARAMETRIC 
X-RADIATION  PRODUCTION 


by 

William  M.  Buckingham 
Lieutenant,  United  States  Navy 
B.S.,  University  of  Washington,  1987 

and 


Authors: 


Approved  by: 


Lisa  R.  Ivey 

Lieutenant,  United  States  Navy 
B.S.,  University  of  Tennessee,  1988 

Submitted  in  partial  fulfillment 
of  the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  PHYSICS 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 


William  M.  Buckingh 


Accesion  For 


NTIS  CRA&I 
OTIC  TAB 
Unannounced 
Justification 


□ 


By . . 

Distribution  / 


Availability  Codes 


Dist 


tL 


Avail  and/or 
Special 


William  B.  Colson,  Chairman,  Department  of  Physics 


iii 


iv 


ABSTRACT 


Previous  measurements  of  parametric  x-radiation  (PXR)  spectra  from  mosaic 
graphite  observed  yields  that  greatly  exceeded  theoretical  predictions.  At  the  Bragg 
scattering  condition,  the  interaction  of  relativistically  charged  particles  and  the  atomic 
planes  of  a  crystal  produces  scattering  of  virtual  photons.  This  is  known  as  PXR.  The 
discrepancies  between  actual  yields  and  theoretical  predictions  might  possibly  have  been 
caused  by  the  mosaic  spread  of  the  atomic  planes  of  the  crystal.  Re-analysis  of  previous 
data  has  shown  a  miscalculation  of  the  effective  target  thickness  used  in  earlier  work. 
There  was  also  an  error  in  the  calculation  of  the  theory  that  was  used  for  comparison. 
Correction  of  the  two  errors  brought  the  data  and  theory  in  closer  agreement,  but  the 
mosaic  spread  does  change  the  intensity  of  the  PXR. 

In  this  thesis,  three  samples  of  thick  mosaic  carbon  graphite  crystal  with  differing 
mosaic  spreads  are  compared  to  determine  the  effect  of  the  mosaicity  on  PXR  yield.  By 
obtaining  “rocking  curves”  where  the  target  orientation  angle  is  changed  with  respect  to 
the  incident  electron  beam  and  observation  angle,  it  was  found  that  the  mosaic  spread  does 
have  an  effect  on  the  PXR  spectral  distribution.  This  effect  is  quite  significant  and  can  be 
seen  when  viewing  the  “rocking  curves”  of  the  samples.  As  the  mosaicity  increased,  the 
angular  distribution  was  spread  out  and  the  coherence  condition  was  relaxed. 
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I.  INTRODUCTION 

Scattering  of  virtual  photons  produced  by  the  interaction  of  a  relativistic  charged 
particle  and  a  crystalline  material  at  or  near  the  Bragg  condition,  2dsin0=nX,  results  in  the 
production  of  parametric  x-radiation  (PXR).  [Ref  1]  This  scattering  is  associated  with 
the  Coulomb  field  of  the  relativistic  particle  and  the  atomic  planes  of  the  crystal;  in  this 
case,  a  thick  mosaic  carbon  graphite  crystal  was  used.  When  the  incident  electron  beam 
angle  with  respect  to  the  crystal  axis  satisfies  the  Bragg  condition,  virtual  photons  are 
placed  on  the  mass  shell  and  scattered  as  real  photons.  The  first  PXR  theory  was 
developed  by  Ter-Mikaelian  (1971)  as  resonant  radiation  that  was  produced  in  a  thin 
crystal.  [Ref  1]  Development  of  the  theory  for  a  thick  crystal  followed  with  experiments 
in  the  mid  1980’s.  [Ref  2]  In  the  early  1990’s  the  first  PXR  experiments  outside  of  the 
former  Soviet  Union  were  performed  at  Naval  Postgraduate  School  (NPS),  Monterey, 
California.  The  first  experiments  using  a  carbon  target  at  NPS  were  reported  in  June  1992 
[Ref  3]  and  in  December  1992.  [Ref  4]  These  experiments  used  a  compression  annealed 
pyrolytic  graphite  (CAPG)  with  a  mosaic  spread  of  0.45°. 

The  recent  theoretical  work  has  been  done  by  Rule  et  al  at  the  Naval  Surface 
Warfare  Center  in  Silver  Spring,  Maryland.  [Ref  5]  This  theory  predicts  experimental 
yields  due  to  the  influence  of  the  mosaic  spread  in  an  attempt  to  explain  the  discrepancies 
seen  between  previous  experiments  and  theories. 

In  this  thesis,  thick  mosaic  graphite  crystals  with  mosaicities  of  0.45°,  1.31°,  and 
2.5°  were  used  with  a  90  MeV  electron  beam  at  a  Bragg  angle  of  22.5°  which  yielded  x- 
rays  in  the  4  to  30  keV  range.  Analysis  of  these  spectra  confirmed  that  mosaic  spread 
does  indeed  have  an  effect  on  the  angular  distribution  and  coherence  conditions. 
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n.  THEORETICAL  BACKGROUND 

PXR  is  an  important  source  of  x-rays  because  of  the  high  energy  photons  that  are 
produced  with  the  relatively  low  electron  energy  requirements  of  the  experiments.  PXR  is 
produced  by  the  interaction  of  relativistically  charged  particles  and  the  atomic  planes  of  a 
crystal.  When  the  Bragg  condition,  2dsin0=nX,  is  satisfied,  virtual  photons  are  placed  on 
the  mass  shell  and  scatter  out  as  real  photons.  Ter-Mikaelian  was  the  first  to  produce  a 
theory  for  PXR  in  a  thin  crystal.  The  first  theory  for  a  thick  crystal  was  by  Feranchuck 
and  Ivanshin.  [Ref  6]  The  emitted  photons  spectral  and  angular  distributions  can  be 
expressed  by: 


d^N  _  g'  I  |2 
dn^da  2it  “ 


\K,,QJV  +  t\ 


-exp(-^)]5(<?), 


(1) 


where  ^  =  +  (2) 

where  N  is  the  number  of  photons  per  electron,  e  is  the  electron  charge,  0s  is  the  multiple 
scattering  angle,  La=[©Im(go)]'\  go  and  g^  are  the  perpendicular  and  parallel  Fourier 
components  of  the  dielectric  susceptibility,  x  is  the  reciprocal  lattice  vector,  n^  is  the 
photon  index  of  refi'action,  ©  is  frequency  of  the  emitted  photon,  v  is  particle  velocity,  Kj^ 
is  the  component  of  the  wave  vector  perpendicular  to  the  crystal  plane.  Feranchuck  and 
Ivanshin  also  found  that  the  number  of  photons  per  electron  per  solid  angle  per  electron 
energy  can  be  found  using  angular  displacement  as  follows: 

d^N  _  1  g'  1  I  ,2  0,'cos^(20g)  +  e/ _ 1 _ 

dd.ddydo)  le  he  Ws  {y-^  +0/  +0/  +\g,\)  i6sin"05(^-|f)'  +p/  ’ 

where  0x  and  0y  represents  the  angular  displacement  from  the  Bragg  condition,  where  ps  is 
defined  as  X/[La(©)27c]  with  La(©)is  the  photon  absorption  length,  cob  is  the  frequency  of 
the  emitted  photon  at  Bragg  condition  ,  y  is  defined  as  l/Vl-(v/c)^,  A©  is  the  difference 
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between  the  emitted  photon  o  and  (Ob-  [Ref.  7]  To  obtain  the  angular  distribution  of  the 
spectral  PXR  spectra,  Equation  (3)  can  be  integrated  over  ©. 

This  theory  when  used  by  DiNova  [Ref  4]  yielded  values  of  N  that  were  much 
smaller  than  the  experimental  photon  production  rate  for  all  spectral  orders  and  for  each 
angle  at  or  near  the  Bragg  condition  for  the  “near  field”  case  and  values  that  were  larger 
thqn  experimental  yields  for  the  far  field  case.  Table  1  shows  a  sampling  of  values  from 
DiNova  and  the  theoretical  values  produced  by  the  above  equation  integrated  over  co. 
There  were  corrections  made  to  the  values  calculated  by  DiNova  in  the  value  of  the 
effective  tin  thickness  as  will  be  presented  in  the  data  analysis  section. 


Table  1.  Comparison  of  Past  Experimental  Data  Corrected  for  Effective  Tin  Thickness 

with  Theoretical  Calculations. 

Rule’s  calculations  presented  at  Tomsk  by  Rule  et  al  [Ref  5]  convolves  PXR 
spectral-angular  distributions  with  Gaussian  models  of  the  mosaic  tilt  angle  distributions 
and  the  effects  of  mosaic  spread,  a™,  is  brought  into  play.  The  effect  of  multiple 
scattering,  Os ,  is  also  included  in  the  following; 
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and  the  effects  of  mosaic  spread,  a^,  is  brought  into  play.  The  effect  of  multiple 
scattering,  CTj  ,  is  also  included  in  the  following: 


with  Cx  =  <5-^  =  Cm  +  used  to  add  the  mean  square  mosaic  tilt  angle  to  the  mean 
square  scattering  angle,  ©xi,  0x2,  0yi  and  ey2  are  the  limits  of  integration  and  in  this  case 
are  all  equal  to  the  radius  of  the  circular  detector  aperture  in  radians,  (DB(n)  =  n7Cc/dsin(0B) 
where  d  is  the  interplanar  spacing,  y  is  the  Lorentz  factor,  Xo  is  the  mean  dielectric 
susceptibility,  and  Xio  is  the  dielectric  structure  factor  This  theory  appears  to  more 
accurately  support  the  experimental  results  obtained. 

The  effects  of  finite  aperture  size  must  be  considered  when  working  with  PXR. 
These  effects  depend  on  the  aperture  area  and  the  distance  between  the  target  and  the 
detector.  From  the  above,  the  solid  angle  of  view  for  the  detector  can  be  determined  and 
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this  can  be  used  to  determine  how  much  of  the  PXR  is  actually  being  viewed  by  the 
detector. 

There  are  two  regions  for  a  finite  detector  aperture  called  the  narrow  line  limit 
regime  and  the  aperture  limited  bandwidth  regime.  [Ref  4]  For  the  first  case,  A9x=Ps,. 
For  the  second  case,  A0x»ps  and  the  distribution  can  be  approximated  as 


(9) 


(10) 


and  tty  =  where  ’ 

Bph  =  m^/E^  +  0s^  + 1  go  1 .  This  is  the  regime  with  which  this  thesis  will  be  concerned. 

Within  this  regime,  the  “near  field”  is  defined  as  A0x  >  FWFIM  of  J2(u)  and  the 
“far  field”  is  defined  by  A0x  <  FWHM  of  JzCu).  The  “near  field”  and  “far  field”  cases  are 
presented  in  Figures  (1)  and  (2),  respectively.  For  this  thesis,  the  experiments  were  done 
in  the  “near  field”,  however,  there  are  some  references  to  previous  “far  field”  expenments 
and  results.  As  can  be  seen  in  Figure  (1),  in  the  “near  field”  there  is  very  little  finite 
aperture  effects  and  most  of  the  PXR  is  within  the  field  of  view  of  the  detector.  On  the 
other  hand,  care  must  be  taken  in  the  analysis  of  the  “far  field”  data  as  a  large  portion  of 
the  PXR  is  out  of  the  view  of  the  detector  due  to  its  limited  solid  angle  size. 
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Figure  1  .  h  and  the  field  of  view  for  the  “near  field”  case 


Figure  2  .  J2  and  the  field  of  view  for  the  “far  field”  case. 
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m.  PXR  EXPERIMENTAL  DESCRIPTION 
The  Naval  Postgraduate  School  Electron  Linear  Accelerator  (LINAC)  was  used  to 
accelerate  dark  current  electrons  to  energies  of  approximately  95  MeV.  [Ref  7]  These 
electrons  then  struck  three  separate  targets  of  compression  annealed  pyrolytic  graphite 
(CAPG)  of  differing  mosaicities.  Initial  energy  calibration  was  performed  using 
sandwiched  foils  of  tin  (Sn)  and  Yttrium  (Y).  In  order  to  provide  in  situ  energy  calibration 
and  the  ability  to  calculate  photon  yields  [Ref  4],  a  tin  (Sn)  foil  was  placed  behind  the 
samples.  “Rocking  curve”  data  was  obtained  by  placing  the  target  so  that  its  crystal  plane 
was  placed  at  an  angle  of  21°  with  respect  to  the  electron  beam  direction  and  then  rotating 
the  target  in  0.5°  increments  from  the  electron  beam  until  26°.  This  allows  for  comparison 
of  the  relative  intensities  of  the  n*  order  peaks  between  samples. 

A.  ELECTRON  BEAM  ALIGNMENT  AND  STEERING 

The  LINAC  is  shown  in  Figure  3.  All  operations  are  performed  and  monitored 
from  the  control  room.  Remote  cameras  are  used  in  the  end  station  to  assist  in  obtaining 
proper  electron  beam  alignment.  The  accelerator  consists  of  three  10  ft.  sections,  each 
with  a  separate  Klystron  connected  by  waveguides.  The  LINAC  is  operated  at  2856  MHz 
with  a  pulse  repetition  rate  of  60  Hz.  The  beam  macro  structure  length  is  approximately  1 
fisec.  Thus,  the  LINAC  duty  cycle  requires  the  ability  to  count  at  60  MHz  when  a  single 
photon  is  detected  during  each  LINAC  pulse.  [Ref  4]  In  order  to  satisfy  the  requirement 
for  low  count  rate,  the  LINAC  was  adjusted  during  operation  to  maintain  one  count  every 
five  to  ten  machine  macropulses.  This  also  ensured  that  simultaneous  counts  of  two 
photons  would  not  appear  as  one  count  of  a  higher  energy  photon.  This  is  of  concern  as 
PXR  peaks  are  integral  multiples  of  each  other  with  respect  to  energy  and  could  result  in 
“pile  up.” 

A  reference  laser  and  retro  optics  is  used  to  position  the  target  ladder  on  which  all 
target  samples  are  mounted.  A  beam  is  then  established  in  the  accelerator  after  a  vacuum 
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is  drawn  in  the  target  chamber.  The  quadrupole  magnets  are  used  to  provide  coarse 
steering  of  the  electron  beam  towards  the  target  Fine  adjustment  of  beam  position  is 


Klystron  Gallery  i 


Figure  3.  Overhead  view  of  the  Naval  Postgraduate  School  100  MeV  Linear  Accelerator 
(LINAC). 

obtained  by  steering  the  beam  through  a  small  pinhole  in  the  target  ladder  using  the  end 
station  steering  magnets 

Target  chamber  geometry  is  shown  in  Figure  4.  The  alignment  laser  is  located  90° 
from  the  entrance  port.  The  Canberra  Si(Li)  detector  [Ref.  8]  is  located  135°  from  the 

^ - --^target  chamber 


incoming  entrance 

1  X 

rport 

exit  port 
_ ^  .1— 

electron  - 7 

beam 

V 

"x.  ) 

1  /  1 

phosphorescent 

screen 

detector 

alignment  laser 

Figure  4,  Overhead  view  of  LINAC  end  station  experimental  setup. 


entrance  port.  The  target  ladder  could  be  raised  and  lowered  remotely  and  it  consisted  of 
a  two  axis  goniometer  to  allow  rotation  of  the  mosaic  crystal.  This  same  detector 
geometry  was  maintained  for  all  data  runs. 

B.  SPECTROSCOPY  ELECTRONICS 

As  with  all  detectors,  the  Canberra  Si(Li)  detector  is  very  susceptible  to  ground 
loops  as  well  as  radiated  noise  from  the  klystrons.  Thus,  the  end  station  is  enclosed  in  a 
metal  mesh  cage  to  reduce  the  interference  of  radiated  noise  from  other  sources.  Lead 
and  borated  polyethylene  bricks  were  also  used  around  the  Si(Li)  detector  to  reduce  the 
background.  Figure  5  shows  the  electronics  setup.  All  system  components  were  located  in 
the  control  room  with  the  exception  of  the  detector  and  amplifier  which  were  located  in 
the  LINAC  end  station. 


_  TENNELEC 

SifLD  detector  ~^TC  244 
Amplifier 


TENNELEC  TC 
■^308  Dual  Linear 
Gate 

- 7K - 


Figure  5.  Electronics  setup. 

The  signal  received  by  the  Si(Li)  detector  is  amplified  by  a  Tennelec  TC  244  amplifier.  A 
Tennelec  TC  308  dual  linear  gate  was  used  to  reduce  unwanted  signal  by  providing  a 
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25  |j.sec  gate  during  which  signal  could  be  received.  This  effectively  made  the  counting 
system  only  active  during  the  LINAC  machine  pulse.  Triggering  of  the  gate  was 
performed  using  a  Stanford  Research  System  Inc.  model  DG  535  digital  delay/pulse 
generator.  By  introducing  a  delay  time  (To)  with  respect  to  the  machine  s  start  sequence 
time  (To),  the  gate  was  adjusted  to  coincide  with  the  arrival  of  the  beam  pulse,  thus,  To 
was  adjusted  to  46  psec. 

Another  problem  inherent  in  energy  calibration  of  the  spectra  is  the  affect  of 
klystron  noise  on  the  signal.  The  LINAC  produces  60  pulses  of  one  microsecond  duration 
every  60  seconds.  This  gating  has  no  effect  on  the  calibration  spectrum  and  assists  in 
limiting  noise  which  is  non-coincident  with  PXR  and  Sn  fluorescence  production.  On  the 
other  hand,  klystron  noise  adds  a  negative  bias  to  the  received  detector  signal  which 
produces  a  negative  bias  to  the  calibration  curve  if  one  relies  solely  on  radioactive  sources 
for  energy  calibration.  Figure  6  shows  the  klystron  noise  signal  in  relation  to  detection 
system  gating  For  this  reason  calibration  spectra  were  taken  prior  to  each  day’s  data  run 
using  Sn  and  Y  to  eliminate  this  potential  source  of  error 


Figure  6.  Representative  Shape  of  Klystron  Noise  During  PXR  Data  Collection.  Line 

Indicates  Approximate  Gated  Portion. 


C.  ENERGY  CALIBRATION 

During  the  initial  energy  calibration  conducted  each  day  prior  to  the  rocking 
curves  the  LINAC  was  operated  at  full  gun  grid  voltage.  The  beam  strikes  sandwiched 
foils  of  tin  (Sn)  and  Yttrium  (Y)  in  order  to  observe  the  and  Kp  fluorescent  lines. 
Energy  of  the  Y  Kx  is  14  957  keV  and  the  energies  of  the  Sn  Kx  and  Kp  lines  are 
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25.270  keV  and  28.483  keV  respectively.  [Ref.  9]  PCA-II  pulse  height  analyzing 
software  [Ref  10]  was  used  to  ensure  calibration  of  all  512  channels.  “Rocking  curve” 
data  was  then  taken  beginning  at  21.0°  from  the  electron  beam  and  proceeding  in  0.5° 
increments  to  26.0°  with  the  LINAC  operating  at  dark  current.  This  is  required  in  order 
to  limit  the  photon  count  rate  since  the  photon  detector,  a  Canberra  model  SI200250 
Si(Li)  detector,  has  a  resolution  time  of  12  psec  determined  by  the  pulse  shaping 
preamplifier.  Dark  current  is  the  term  used  to  describe  operation  of  the  LINAC  with  no 
gun  grid  voltage  applied.  Thus,  only  stray  electrons  are  accelerated  toward  the  target 
crystal. 

In  order  to  provide  in  situ  energy  calibration,  a  tin  (Sn)  foil  was  placed  behind  each 
CAPG  sample.  Table  2  provides  target  and  tin  (Sn)  thickness  for  each  sample.  Each 
thickness  measurement  was  obtained  by  taking  the  average  of  five  thickness 
measurements  at  different  positions  on  the  sample.  Each  target  crystal  has  a  density  (p)  of 
2.260  ±  0.005  grams  per  cubic  centimeter.  The  {002}  reflecting  plane  spacing  is  3.356  ± 
0.003  Angstroms. 

The  Union  Carbide  specifications  for  the  mosaicities  of  the  samples  were  0.4±0.1°, 
1. 2+0.2°,  and  3. 5+1. 5°.  Since  the  tolerances  from  the  prescribed  mosaic  spread  were 
quite  large,  the  samples  were  sent  to  Naval  Surface  Warfare  Center  (NSWC)  for 
determination  of  the  actual  mosaic  spreads.  This  was  measured  with  an  x-ray 
diffractometer.  The  resulting  spectra  from  the  x-ray  diffraction  are  presented  in  Figures  7 
through  9.  [Ref  1 1]  For  each  sample,  the  mosaic  spread  was  determined  by  the  full  width 
at  half  maximum  (in  degrees)  of  these  curves.  The  resulting  spreads  were  0.45°,  1.2°, and 
2.5°. 


Mosaic  Spread 

0.45° 

1.31° 

2.5° 

Target  Thickness 

1389  +  1  pm 

1659  ±  1  pm 

1789+  1  pm 

Sn  Thickness 

152  ±1  pm 

1 54  ±  1  pm 

156  ±  1  pm 

Table  2.  Sample  and  Sn  Thickness  for  each  Sample. 
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_  Energies  (in  keV)  for  each  channel  were  established  during  the  calibration  runs  by 
the  PCA-II  pulse  height  analyzing  software  using  the  following  equation: 

E  =  ao  +  ai(channel)  +  a2(channel)^  ( 12) 

where  ao,  ai  and  aa  are  constants  determined  by  performing  Gaussian  fits  of  the  counts 
received  in  each  region  of  interest  surrounding  the  Yt  and  Sn  calibration  peaks.  These 
constants  for  each  series  of  data  is  shown  in  Table  3.  The  energy  calibration  was  verified 
by  comparing  the  and  lines  produced  by  the  Sn  backing  used  during  each  PXR  data 
run  against  the  known  value.  Good  agreement  was  observed  for  all  PXR  spectral  data. 


Mosaic  Spread 

0.45® 

1.31® 

2.5® 

ao 

1.95 

1.86 

1.34 

ai 

0.0722 

0.0712 

0.0764 

32 

0 

0 

8.53e-6 

Table  3.  Calibration  Constants  for  each  Sample. 
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Counts 


30000 


Rocking  Curve  Mosaic  Graphite  ZYA 

Copper  Ka  ( 8  keV ) 


Figure  7.  Mosaicity  Measurement  curve  for  0.45°  Sample. 
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Counts 


Rocking  Curve  Mosaic  Graphite  ZYD 

Copper  Ka  ( 8  keV ) 


Figure  8.  Mosaicity  Measurement  curve  for  1.31°  Sample. 
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IV.  DATA  AND  ANALYSIS 


A.  CALCULATION  OF  PXR  PEAK  AREAS 

The  data  from  the  detector  was  collected  using  a  PCA-II  pulse  height  analyzer 
[Ref  10]  in  ASCII  format.  Raw  data  spectra  for  the  0.45°.  1.31°.  and  2.5°  mosaic 
samples  are  shown  in  Figures  10.  11.  and  12  respectively.  Data  was  recorded  for  each 
sample  spectrum  for  approximately  30  minutes  with  electron  energies  (Ee)  of  99  05  MeV, 
95.32  MeV,  and  97.62  MeV  respectively.  The  spectra  were  then  input  into  Quattro  Pro 
(QPRO)  spreadsheets  [Ref  12]  so  that  it  would  be  in  proper  format  for  the  spectrum 
analyzing  software.  Peakfit  v3.0  [Ref  13]  was  used  to  remove  the  background  noise  and 
produce  a  Gaussian  fit  for  each  spectral  peak.  Figure  13  a)  b)  and  c)  show  the 
uncorrected,  background  removed,  and  attenuation  corrected  spectrums  respectively  for 
the  0  45°  sample  at  22.5°.  The  program  also  integrated  each  peak  to  produce  a  net  area 
for  each  peak  and  the  full  width  at  half  maximum  values  for  each  peak  of  the  spectrum. 
The  values  were  placed  in  a  Nficrosoft  Excel  worksheet  [Ref  14]  and  corrected  for 
attenuation  using  Photcoef  [Ref  15]  Photcoef  was  used  to  calculate  the  attenuation 
coefiScients  for  the  Kapton  (C22O5N2H10)  window,  the  Berylium  (Be)  window  and  the  air 
gap  at  each  energy.  Thicknesses  for  each  of  these  were  0  0025  cm  0.005  cm,  and  0.5  cm 
respectively.  At  low  energies,  the  attenuation  coefficients  varied  across  each  peak.  For 
the  first  peaks  only,  the  attenuation  coeflSdents  for  the  left  most  channel,  the  peak  center 
and  the  right  most  channel  were  added  in  quadrature  and  averaged  to  obtain  a  more 
accurate  correction  for  the  peak.  For  the  higher  order  peaks,  the  attenuation  coefficients 
were  approximately  constant  across  the  peak  and  the  coefficient  was  obtained  for  the 
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center  channel  of  the  peak  and  applied  to  the  entire  peak.  The  areas  were  then  corrected 
to  account  for  these  attenuations  as  follows: 


(13) 


where  No  is  the  total  number  of  photons  before  attenuation.  N  is  the  total  number  of 
photons  after  attenuation,  p  is  the  attenuation  coefficient,  p  is  the  density  of  the 
attenuating  medium,  and  t  is  the  distance  travelled  in  the  attenuating  medium.  The 
efficiency  of  the  detector  had  to  also  be  accounted  for  in  the  corrections.  The  efficiency 
curve  for  the  Canberra  detector  at  various  energies  can  be  found  in  Figure  14.  [Ref  8] 
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Figure  10.  Raw  0.45°  Mosaicity  Spectra  at  Ee  =  99.05  MeV 
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Figure  11.  Raw  1.3 1°  Mosaicity  Spectra  at  =  95.32  MeV 
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Figure  12  Raw  2  5°  Mosaicity  Spectra  at  =  97.62  MeV. 


25 


Energy  (keV) 

Figure  13.  a)  Uncorrected  PXR  spectrum  for  0.45°  mosaic  spread  at  22.5°.  b)  The  same 
spectrum  with  background  removed,  c)  Spectrum  corrected  for  attenuation  between  the 
sample  and  the  detector. 
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Figure  14.  Efficiency  Curves  for  Canberra  Si(Li)  Detector. 

B.  CALCULATION  OF  PXR  ABSOLUTE  YIELD 

The  corrected  PXR  areas  are  then  used  as  follows  to  calculate  the  number  of 


photons  produced  for  each  electron. 


Combining  Equations  14  and  15  yields: 


Area^xR^  TIH  pNJiJ 

(^Arecij^  ^A 
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where  AreapxR  and  Area-nN-  are  the  integrated  intensities  of  the  respective  peaks,  corrected 
for  both  attenuation  and  detector  efficiency,  AWtin  is  the  atomic  weight  of  tin  a  is  the 
electron  interaction  cross  section  (1  Oe-26  m^)  [Ref  16]  Qdet  is  the  detector  solid  angle 
(4.15e-4  sr.),  fdex  is  the  de-excitation  transition  probability  (0.712)  [Ref  17],  and  t  is  the 
effective  tin  thickness.  The  effective  thickness  used  in  the  above  equation  was  calculated 
by  dividing  the  thickness  of  the  tin  (ti)  by  the  sine  of  the  angle  of  incidence  of  the  electron 
beam.  The  thickness  of  the  tin  (ti)  was  152  ±  1  pm,  154  ±  1  pm,  and  156  +  1  pm  for  the 
0.45°,  131°,  and  2.5°  mosaic  samples  respectively.  The  geometry  of  the  electron  beam  as 
it  strikes  the  Sn  foil  presented  in  Figure  15.  This  calculation  was  done  for  the  first  peak  at 
each  angle  of  the  rocking  curve  for  all  three  samples  The  number  of  photons  per  electron 
(N)  versus  angle  for  the  first  and  second  order  PXR  peaks  are  presented  in  Figures  16  and 
17  for  each  sample. 


electron  beam 


Figure  15.  Configuration  in  Target  Chamber  Used  to  Determine  Tin  Effective  Thickness, 
ti  is  the  Normal  Target  Thickness.  0b  is  the  Bragg  Angle  or  the  Angle  Between  the 
Target  Face  and  the  Electron  Beam  Direction. 


29 


Figure  16,  Comparison  of  Corrected  First  Order  PXR  Intensities  (N).  0.45°,  1.31°,  and  2.5°  Refer  to  the  Mosaic  Spread  of  the 
Graphite  Targets.  The  Electron  Energies  (Ee)  are  99.05  MeV,  95.32  MeV,  and  97.62  MeV  Respectively. 


N(photons/electron)x  10''-7 
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SECOND  PEAK  COMPARISON 
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C.  CALCULATION  OF  INTENSITY  RATIOS 


The  higher  order  peaks  were  compared  to  the  first  peak  absolute  5deld  for  each 
angle  of  the  samples  to  give  an  intensity  ratio.  These  calculations  are  presented  in  Table 
4  and  will  be  used  for  analysis  of  theoretical  predictions.  For  each  sample,  the  values 
presented  in  Figure  16  were  normalized  to  the  highest  first  order  intensity  to  provide  a 
comparison  of  the  relative  intensities  of  the  first  order  PXR  peaks  across  the  “rocking 
curve.”  These  values  are  presented  in  Figure  18.  The  same  comparison  was  made  for  the 
second  order  PXR  peaks  and  is  presented  in  Figure  19. 

In  order  to  provide  a  comparison  of  the  data  with  theory  for  all  three  samples, 
Equation  3  was  integrated  over  q,  0x,  and  0y  for  the  solid  angle  subtended  by  the  detector. 
These  calculations  are  referred  to  here  as  theory  without  scattering.  To  accovmt  for 
multiple  scattering  and  mosaic  spread.  Equation  17  was  used.  [Ref  4]  These  calculations 


ci‘N  ^  f  r,  _  I'llkla’jf  + 

3eje,  S4;r  c  ‘ sin'-e,  + 

2 

where  0^/.^  =^+0/ +|^o| 

and  0,~  —  nwiacai  0^ beam-dn  “t"  0^ mos 


(17) 

(18) 
(19) 


are  referred  to  as  theory  with  scattering.  Appendbc  A  shows  a  sample  of  the  MathCad 
[Ref  18]  program  used  to  calculate  the  values  shown  for  theoretical  calculations,  with  and 
without  taking  into  account  the  effects  of  mosaic  spread  and  multiple  scattering.  Figures 
20  through  25  give  a  comparison  of  the  data  with  the  resulting  theoretical  curves  for  each 
sample,  including  both  the  first  and  second  order  PXR  peaks.  Agreement  between  the 
and  theory  appears  to  break  down  as  the  mosaic  spread  increases.  Thus,  the  method 
presented  by  Rule  et  al  of  convolving  the  PXR  spectral-angular  distributions  with 
Gaussian  models  of  the  mosaic  tilt  angle  appears  to  provide  better  agreement.  Table  5 
provides  the  theoretical  intensities  and  ratios  for  the  first  and  second  order  PXR  peaks  at 
Bragg  condition  provided  by  Rule  for  a  detector  solid  angle  of  4.15e-4  sr.  [Ref  19] 
Further  calculations  using  this  model  are  in  progress. 
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Figure  20.  Corrected  First  Order  PXR  Intensity  Comparison  with  Theory  for  0.45°  Mosaic  Sample. 


N(photons/electron) 
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0.45  DEG.  FIRST  PEAK  DATA  VS.  THEORY 
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Figure  21 .  Corrected  Second  Order  PXR  Intensity  Comparison  with  Theory  for  0.45°  Mosaic  Sample. 


Figure  22.  Corrected  First  Order  PXR  Intensity  Comparison  with  Theory  for  1.31°  Mosaic  Sample. 


N(photons/electron) 


38 


131  DEG.  FIRST  PEAK  DATA  VS.  THEORY 


1.31  DEG.  SECOND  PEAK  DATA  VS.  THEORY 
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Figure  23.  Corrected  Second  Order  PXR  Intensity  Comparison  with  Theory  for  1.31°  Mosaic  Sample. 
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2.5  DEG.  FIRST  PEAK  DATA  VERSUS  THEORY 


2.5  DEG.  SECOND  PEAK  DATA  VS.  THEORY 


O) 

c 

c 

£ 

(0 

i 

o 

Q 

c/3 

O 

o 

C/3 

> 

> 

>. 

o 

O 

o 

x: 

£ 

1- 

(U0J}39i3/SU0)0(4d)N 


41 


Figure  25.  Corrected  Second  Order  PXR  Intensity  Comparison  with  Theory  for  2.5°  Mosaic  Sample. 


Mosaic  Spread 

II 

N  =  2 

0.45° 

2.14e-5 

7.81e-6 

0.36 

1.31° 

1.63e-5 

5.79e-6 

0.36 

2.5° 

1.06e-5 

3.70e-6 

0.35 

Table  5.  Theoretical  Intensity  Comparisons  Using  Rule  et  al  for  Detector  Solid  Angle  of 

4.15e-4  sr.  [Ref  19] 

D.  CORRECTIONS  TO  PREVIOUS  DATA 


The  same  0.45°  mosaic  spread  sample  was  used  in  recent  experiments  as  was  used 
in  prior  experiments.  DiNova  [Ref  4]  measured  photon  yield  (N)  with  the  detector  placed 
at  100  cm  and  29  cm  from  the  target  which  will  be  referred  to  as  the  “far  field”  and  “near 
field”  cases  respectively.  DiNova  calculates  the  effective  thickness  of  tin  used  in  the 
absolute  yield  calculations  by  dividing  the  thickness  of  the  tin  sample  by  the  cosine  of  0b. 
[Ref  16]  As  presented  earlier,  this  should  be  divided  by  the  sine  of  0b.  Once  this 
correction  is  made  to  DiNova’s  data,  the  comparison  between  her  data,  and  recent  data 
reveals  similar  results.  Differing  detector  solid  angles  account  for  most  of  the  difference 
between  recent  data  and  the  older  near  and  far  field  data.  Other  differences  include 
varying  tin  foil  thicknesses,  different  detectors,  and  different  data  analysis  procedures. 
Table  6  summarizes  these  differences.  In  order  to  determine  how  large  the  effect  of  the 
differing  data  analysis  techniques  was,  DiNova’s  raw  data  was  reprocessed  in  the  exact 
manner  as  for  the  three  samples  presented  here.  For  comparison  purposes.  Table  7  shows 
the  first  order  peak  data  using  the  corrected  effective  tin  thickness  from  DiNova’s  thesis 
and  DiNova’s  data  re-analyzed.  Figures  26  and  27  compare  this  data  graphically.  An 
attempt  was  made  to  re-analyze  the  data  for  the  25.6°  spectrum,  but  the  first  order  peak 
was  not  distinguishable.  DiNova’s  re-analyzed  near  field  and  far  field  and  the  recent  0.45° 
mosaic  sample  data  are  also  con;  ired  in  Table  8  to  show  how  the  distance  from  detector 
to  target  effects  the  absolute  yield.  It  is  important  to  note  that  the  new  data  fell  at  a 
distance  that  was  in  between  DiNova’s  far  and  near  field  results. 
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Near  Field  Data 

New  Data 

Far  Field  Data 

2.39e-3  sr. 

4.15e-4  sr. 

2.0e-4  sr. 

Tin  Thickness  (ti) 

Detector 

Analysis  Procedure 

Summed  Counts  in 
PXR  Peak 

Peakfit  v3.0 

Summed  Counts  in 
PXR  Peak 

Table  6.  Summary  of  Differences  between  Prior  and  Recent  Experiments. 


Angle  (0b) 

Near  Field 
Data 

Near  Field  Data 
Re-Analyzed 

Far  Field 
Data 

liiiM 

20.1° 

3.80e-6 

3.60e-6 

- 

- 

21.1° 

1.36e-5 

1.28e-5 

2.27e-6 

5.12e-7 

21.6° 

3.01e-5 

2.73e-5 

5.46e-7 

9.30e-7 

22.1° 

4.70e-5 

2.71e-5 

1.14e-6 

2.55e-6 

22.6° 

6.22e-5 

5.49e-5 

4.68e-6 

8.49e-6 

23.1° 

3.82e-5 

3.30e-5 

3.87e-6 

2.07e-6 

23.6° 

1.58e-5 

1.28e-5 

8.74e-7 

1.59e-6 

24.1° 

8.07e-6 

7.33e-6 

4.25e-7 

7.75e-7 

24.6° 

6.47e-6 

4.77e-6 

3.06e-7 

4.59e-7 

25.6° 

4.76e-6 

- 

- 

- 

Table  7.  Comparison  of  DtNova’s 

Data  Corrected  for  ti 

le  Effective  Tin  Thickness  and 

DiNova' s  Data  Re- Analyzed  Using  PeakFit  (0.45°  Mosaic  Spread) 


By  changing  the  distance  from  the  target  to  the  detector,  the  solid  angle  of  the 
detector  is  changed.  As  can  be  seen  in  Figure  28  for  the  three  different  solid  angles,  there 
is  a  marked  effect  on  the  photon  5aeld  (N)  as  the  distance  is  increased 
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1 

Near  Field  Data 

New  Data 

Far  Field  Data 

2.39e-3  sr. 

4.15e-4  sr. 

2.0e-4  sr. 

N(photons/elec.) 

N(photons/elec.) 

N(photons/elec.) 

3.60e-6 

1 

b 

0 

1.21e-6 

-15° 

1.28e-5 

1.50e-6 

5.12e-7 

0 

o 

2.73e-5 

3.29e-6 

9.30e-7 

-0  5° 

2.71e-5 

7.81e-6 

2.55e-6 

0° 

5.49e-5 

2.27e-5 

8.49e-6 

0.5° 

3.29e-5 

1.63e-5 

2.07e-6 

0 

O 

1.28e-5 

9.12e-6 

1.56e-6 

1.5° 

7.33e-6 

3.99e-6 

7.75e-7 

2.0° 

4.76e-6 

2.85e-6 

4.59e-7 

2.5° 

1.56e-6 

3.0° 

1.18e-6 

Table  8.  Comparison  of  0.45°  Mosaic  Sample  Results  by  Solid  Angle. 


As  can  be  easily  seen  from  Figures  16  and  17.  as  the  mosaic  spread  of  the  crystal 
increases,  the  angular  distribution  of  the  “rocking  curve”  also  increases.  The  PeakFit 
software  was  again  used  to  obtain  the  full  width  at  half  maximum  for  the  first  order  peaks. 
The  resulting  distribution  appears  linear  and  is  shown  in  Figure  29  Thus,  increasing 
crystal  mosaicity  relaxes  the  coherence  condition  and  spreads  out  the  resulting  rocking 
curve”  angular  distribution. 
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Figure  27,  Corrected  and  Re- Analyzed  Data  for  Far  Field  0.45°  Mosaic  Sample  for  Detector  Solid  Angle  of  2  Oe-4  sr. 


N(photons/electron) 
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FAR  FIELD  COMPARISON 


SOLID  ANGLE  EFFECTS  ON  PHOTON  YIELD 


^  vO  |.x(uoi)33|3/suo)oqd)N 
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e  28  0  45°  Mosaic  Sample  Intensities  by  Solid  Angle. 


ANGULAR  DISTRI BUTION 


I  I  I  I  I  I  I  I  I  1  I  I  I  I  I- 1— <■ 
0.5  1  1.5  2 

Mosaic  Spread  {deg 


♦  3.6 


I  -{■  -<-  'l  I  I  I 
2.5  3 


Figure  29.  Full  Width  at  Half  Maximum  of  First  Order  Peaks  as  a  Function  of  Mosaic 
Spread. 

E  ERROR  ANALYSIS 

There  are  two  possible  sources  of  errors  in  this  experiment,  namely,  systematic  and 
random  errors.  The  sources  of  systematic  errors  are  caused  by  limitations  of  the  systems 
used  and  effect  all  measurements  in  the  same  fashion.  These  errors  are  often  hard  to 
analyze  and  to  even  detect  by  statistical  analysis  Random  errors  are  experimental 
uncertainties  that  can  be  revealed  by  doing  the  same  measurements  a  number  of  times  to 
get  an  average  value  of  the  measurement.  Treatment  of  the  two  types  of  errors  are 
handled  differently.  [Ref  20] 

In  handling  systematic  errors,  the  best  analysis  is  to  predict  the  possible  sources 
and  verify  that  these  errors  will  be  less  than  the  precision  that  is  required  for  the 
experiment.  The  systematic  errors  for  this  experiment  are  the  angle  of  target  incidence 
which  is  due  to  the  mechanical  motor  that  turns  the  target  and  is  subject  to  over  and/or 
under  rotation.  The  beam  alignment  is  also  a  systematic  error  that  would  not  have  an 
effect  on  each  individual  run,  but  when  comparing  different  days  runs,  this  error  will  be 
evident.  This  error  also  effects  the  number  of  electrons  that  actually  strike  the  target 
crystal  and  tin  foil  which  is  a  value  that  is  calculated  from  the  areas  of  the  spectra.  The 
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combination  of  these  errors  produced  a  shifting  of  maximum  intensities  of  0.5°  for  the 
0.45°  mosaic  sample,  and  1.0°  for  the  1.31°  and  2.5°  samples.  This  error  is  evident  when 
comparing  the  maximum  experimental  intensities  to  the  calculated  theoretical  peak 
intensities. 

The  random  errors  are  handled  by  adding  in  quadrature  the  errors  involved  with 
the  variables  used  in  calculating  the  number  of  photons  per  electron  as  seen  in  Equation 
16.  The  only  two  terms  which  contribute  to  this  calculation  are  the  areas  of  the  PXR 
peaks  and  the  area  of  the  tin  peak.  Equation  20  was  used  to  calculate  the  fractional  error 
in  the  number  of  photons  per  electron.  The  random  errors  in  area  were  analyzed  using 
Equation  21. 


where  F  is  the  counts  in  a  spectral  peak,  Ho  is  the  height  of  the  PXR  peak,  ai  is  the  center 
channel  of  the  PXR  peak,  a2  is  the  standard  deviation  of  the  curve  fitted  peak  used  by 
Peakfit  to  determine  the  best  fit  peak,  and  dao,  dai,  and  da2  are  the  errors  in  each  of  these 
values  respectively  that  are  provided  by  the  Peakfit  Peak  Analysis  Program.  The  integral 
was  processed  using  a  MathCad  [Ref  18]  program  that  is  listed  in  Appendk  B  for  each 
peak  of  the  PX[R  spectrum  and  the  tin  peak  for  each  spectrum.  The  errors  in  the  PXR 
peak  and  the  tin  peak  are  then  added  in  quadrature  to  obtain  the  total  error  for  each 
calculation  of  the  number  of  photons  per  electron.  Table  9  lists  the  results  of  Equations 
20  and  21  as  well  as  those  from  adding  in  quadrature  the  errors  in  the  AreapxR  and  AreanN 
using  Equation  22  for  each  of  the  first  peaks  for  each  spectra.  Figures  30  through  34 

(22) 
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1.45°  Mosaicity _ 1.31°  Mosaicity _ 2.5°  Mosaicii 


MCAD 

iWArea 

MCAD 

1/VArea 

MCAD 

l/vArea 

21.0° 

.0079 

.062 

.0078 

.054 

.0164 

.042 

21.5° 

.0082 

.064 

.0091 

.036 

.0192 

.042 

22.0° 

.0113 

.041 

.0107 

.019 

.0247 

.027 

22.5° 

.0262 

.026 

.0252 

.012 

.0515 

.028 

23.0° 

.1747 

.018 

.0499 

.010 

.0405 

.023 

23.5° 

.0805 

.018 

.0782 

.014 

.0411 

.023 

24.0° 

.0303 

.023 

.0420 

.015 

.0345 

.022 

24.5° 

.0131 

.035 

.0131 

.027 

.0276 

.027 

25.0° 

.0094 

.035 

.0174 

.036 

.0214 

.033 

25.5° 

.0070 

.052 

.0089 

.057 

.0149 

.037 

26.0° 

.0061 

.059 

.0069 

.062 

.0154 

.039 

Table  9.  Estimation  of  the  Random  Error  in  PXR  Intensity  (N)  Calculation  Using 
Equations  20,  21,  and  22. 


show  the  random  error  based  on  Equations  20  and  2 1  in  the  earlier  near  and  far  field  re¬ 
analyzed  0.45°  mosaic  spread  data  as  well  as  the  new  data  for  each  of  the  first  order 
peaks.  Using  this  method,  larger  intensity  peaks  incurred  a  larger  fractional  error  due  to 
the  relatively  larger  size  of  the  peak. 
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0.4S  DEG.  NEAR  FIELD  ERROR 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

gO 

(O 

iO 

e5 

CN 

ivO  LX(UOJp3|3/SUO)OUCl)ty{ 


U) 

o 


o 


O) 

4) 


Bi 

e 


< 


e 

9> 


5 


w 

u 

« 


w 

h 


lO 

d 

I 


LO 


Cs 

cn 

<N 


an 

3 


n 

C/5 

u 

P 

O 

I 


-2 

g 


m 

2 

3 

op 


51 


0,45  DEG,  FAR  FIELD  ERROR 


zs 


i  vO  l■x(uoJu^!^/suoloud)^J 


Figure  31.  Random  Error  Distribution  for  0.45°  Mosaic  Sample  With  Detector  Solid  Angle  of  2.0e-4  sr. 


0t.4.s  degree;  mosaicity 
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1,311  DEGREE  MOSAICITY 


PS 


4-yOIx(oojp3ia/snoioqi)M 


Figure  33.  Random  lirror  Distribution  for  1.31®  Mosaic  Sample  With  Detector  Solid  ^Vn,gle  of  4.1:5e"4  sr. 
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Fiigiure  34,  Random  Error  Distribution  for  2.5°  Mosaicity  Sample  With  Detector  Solid  Angle  of  4. 1 5e-4  5ir. 
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V,  CONCLUSIONS  AND  RECOMMENDATIONS 

As  suspected  from  previous  experimental  and  theoretical  comparisons,  the  mosaic 
spread  of  the  crystal  used  does  have  a  definite  effect  on  the  outcome  of  the  PXR  spectra 
and  the  resulting  intensities  and  absolute  yields.  Much  work  is  being  done  by  Rule  et  al  to 
model  theoretically  the  effect  of  the  mosaic  spread  on  PXR  scattering.  The  theory 
presented  by  Rule  et  al  [Ref  5]  seems  to  predict  within  experimental  errors  the  data 
obtained  through  the  use  of  the  three  samples  for  this  thesis  more  than  the  older  theories 
that  were  previously  used.  The  work  of  Rule  et  al  incorporates  the  mosaic  spread,  which 
the  previous  work  by  Feranchuck  and  Ivansin  [Ref  6]  neglected.  As  can  be  seen  by 
Figures  16,  17,  and  29  an  increase  in  the  mosaic  spread  relaxes  the  coherence  condition 
and  spreads  out  the  angular  distribution  of  the  scattered  photons.  Previous  data  used  a 
miscalculated  effective  thickness  for  the  tin  which  made  the  absolute  calculated  yield 
appear  higher  than  it  should  have  been.  In  the  same  comparison,  the  theory  was  also 
miscalculated  making  it  to  appear  lower.  By  correcting  for  these  two  mistakes,  the  data 
and  the  theory  are  more  comparable,  but  there  still  needs  to  be  the  correction  for  the 
mosaic  spread  as  it  does  have  an  effect  Theoretical  calculations  need  to  be  done  to 
support  the  off  Bragg  condition  measurements  as  well  as  the  higher  order  peak  yields  from 
theory  presented  by  Rule  et  al.  Through  more  experimental  data  and  more  theoretical 
calculations  to  compare  with,  a  better  understanding  of  the  mosaic  effects  on  PXR  spectra 
can  be  achieved. 

Future  thesis  topics  should  investigate  the  same  samples  in  a  “far  field” 
configuration  as  well  as  new  samples  with  different  mosaic  spreads.  Contacts  should  be 
maintained  with  Rule  et  al  to  continue  to  compare  past  and  present  data  with  new 
theoretical  calculations  as  they  become  available. 
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APPENDIX  A.  SAMPLE  THEORETICAL  YIELD  PROGRAM 


Differential  cross  section  calculation 
Carbon  (002) 

Input  parameters: 

degree 
Angstrom 
gm  molar  mass 


Sg  22.5 


d  :=  6.7078 


12.011 


V  :=  35.189 


Za-6 


cubic  angstrom 
Atomic  number 


carbons 

.mcd) 

A- 0.45 

degree 

:=  1659 

p.m 

a  "2-26 

gm/cc 

e  =  95.32 

MeV 

mos 


1.31-7C 

180 


Constants 


:=  1.6021917-10 


19 


Electron  charge 


m:=  9.109558-10' 

c  :=  299792458 


.-31 


a:=  7.297351-10' 


,-3 


Electron  mass 

speed  of  light 

Fine  structure  constant 


‘bar  '■ 


6.626196-10 

(2-n) 


-34 


Plank'S  constant 


N^:=  6.022169-10 


.23 


Avogadro  constant 


Fg  ;=  2.817938-10' 


15 


Classical  electron  radius 


59 


Unit  conversion 


e 


9b 


B 


L:  = 


LIO 


180 

,-6 


sm(eB) 


e  A  =  7t- 


V  :=  v  io 


180 


30 


Pa:=P„-10^  d:=d.l0-^° 


Bragg  frequency  calculation 

i  ,=  4  Order  of  Bragg  reflection 


1-1IC 


(0 


B 


d- sin  (0  b) 


Bragg  frequency  function 


Absorption  length  calculation 
Absptn  :=  READPRN(carbon) 


i  :=  0.. rows(Absptn)  -  1 


Input  is  associated  with  absorption 
data  file  *.prn  generated  by  xcom.exe 
in  units  of  MeV  vs  barns 


Absptnj  Q  :=  log^Abspta;  0)  +  6 


Absptoj  j  :=  log^Abspto;  j)  -  28 


La(®)  •= 


m 


a 


•10 


<o> 


-linterp  Absptn  ,  Absptn 


<1> 


^A-Pa 


Atomic  form  factor  function 

cf.  International  Tables  for  x-ray  crystallography,  Vol.  IV  1974  p71 


II 

■=  1.02  33  :=  1.5886 

a4  :=  0.865 

bj  :=  20.8439 

bj  :=  10.2075  bj  :=  0.5687 

64  :=  51.6512 

0  0=  0.2156 

i:=  1..4 

Input  parameters 

f(x)  :=^aj-exp  -bj-(x-10 


10 


+  c, 


X  should  be  sin(9)/?i 


60 


Dielectric  function 


= 


rg-(2-7Cc)^ 

n-Y-a^ 


sinfeg) 

X(C0)  := - ^ ^ 

2-k-c 

S£-:=16  Input  parameter:  Selection  factor 


Zo(co):=  y(28.8162)- 


■f(x(®)) 


N  calculation 
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56  j  =  0,011522 


50^  :=  0.03-tan ^0  g) 


50^ 


50  2  =  0.011522 


0.511 

Ee 


0  :=  .008727  Z2 


6  :=  .017453  Z3 


0  :=  .02618  Z4 


0  :=  .034907 

Z5 


0  :=  .043633  Zg 


6  =.05236  Z7 


0  :=  .061087  Zg 


1..8 


•6-h50  1 

•50  2 

6-501 

-502 

■  0  -t-  50  1 

•502 

0-501  . 

-502 

•6-h50  1 

•50  2 

_  0  -  50  1 

-50  2 

V0-I-50  1 

•50  2 

0-501  . 

-50  2 

•0-t-50  1 

•50  2 

6-501  . 

-50  2 

•0-1-56  1 

•50  2 

0-  50  . 

J  1  • 

-502 

•0-^501 

•502 

,0-501 

-502 

•0-h50  1 

•50  2 

,  0  -  50  1 

-50  2 

Nb(e^.ev„)de^d0y 

Nb(e^’%)d0^d0^ 

Nb(e^.e^)de^d0^ 

Nb(®^’®v|/)d0^d6y 

Nb(e^-%)d0  4d0^ 
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0.157-Za-(Za  +  l)-.511^-L 


(|x  io(“b)|) 


aPB 


0^^-(cos(2-0B))^  +  e^^ 


0.5n\‘‘  „  2  „  2  „  /-  \2  I  / 

-—  +6^  4-0^  ■"®msqr‘^  (®mos)  "-l^oi^B 


i:=  1..8 


0  :=  .008727 


0  :=  .017453 


0  :=  .02618 


'0  +  50^  r50  2 

0-50^  ^-502 

’0  +  50^  r50  2 

J0-601  ^-802 
(•0  +  501  V50  2 


Nbl(eq.%)de4d0^ 


Nm  0 


0-501  J-502 

■0  +  50  1  r50  2 


Nbi 


0  :=  .034907 


0  :=  .043633 


0-501  J-502 

■0  +  50  1  r50  2 


Nbiie^ 


.0v),)de^d0y 


0-501  J-502 
•0  +  50  1  r50  2 


10-50  1  J-50  2 


Nbl  0 


^,0^)  d0^d0^ 
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0  :=  .05236 


zz. 


e  :=  .061087 

zz. 


i:=4 


z  = 


z  is  w/o 


•0  +  50 

0-  50  1 
•0  +•  50  1 

0-  50  2 


•50  2 

,-50  2 

•50  2 
-50  2 


0 


mos 


1.31-JC 

180 


50^:=  0.011522 
Eg:=  95.32 
L  :=  1659 


zz  = 


scattering,  zz  is  with  scattering 
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APPENDIX  B.  SAMPLE  ERROR  ANALYSIS  PROGRAM 


inputs  for  peak  height;  F  is  the  area  under  the  PXR  peak 

ag  ;=  349.6233  da^,  ;=  3.7878  F  :=  1500.9662 

inputs  for  peak  center  channel; 
aj:=  37.46832  daj:=  .021196 

inputs  for  peak  width; 
a^  :=  1.712699  da^  :=  .02188 
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APPENDIX  C.  DATA  VERSUS  THEORY  FOR  ALL  SAMPLES 
This  appendix  provides  a  complete  comparison  of  the  actual  and  theoretical 
intensity  ratios  between  PXR  peaks.  Theoretical  calculations  without  scattering  are  based 
on  the  integration  of  Equation  3 .  Theoretical  calculations  with  scattering  account  for 
multiple  scattering  and  mosaic  spread  based  on  the  integration  of  Equation  17  Here  we 
also  provide  a  comparison  of  the  observed  intensities  versus  the  theoretical  intensities  both 
with  and  without  accounting  for  multiple  scattering  and  mosaic  spread  as  calculated  using 
the  MathCad  program  found  in  Appendix  A. 
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